On the pinning strategy of complex networks 
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In pinning control of complex networks, a tacit believing is that the system dynamics will be better controlled 
by pinning the large-degree nodes than the small-degree ones. Here, by changing the number of pinned nodes, 
we find that, when a significant fraction of the network nodes are pinned, pinning the small-degree nodes could 
generally have a higher performance than pinning the large-degree nodes. We demonstrate this interesting 
phenomenon on a variety of complex networks, and analyze the underlying mechanisms by the model of star 
networks. By changing the network properties, we also find that, comparing to densely connected homogeneous 
networks, the advantage of the small-degree pinning strategy is more distinct in sparsely connected heterogenous 
networks. 

PACS numbers: 05.45.XL89.75.Hc 



Pinning control refers to the controlling of system dynam- 
ics to a target state by perturbing only partial of the system 
variables Qly]. Due to its convenient operation and high per- 
formance, this technique has been widely used in controlling 
spatiotemporal dynamics in complex systems, including op- 
tical systems JH, plasma 0], neural systems 0], and chem- 
ical and flow turbulence in various systems 0-0]. In pin- 
ning control, one of the focusing issues is about how to op- 
timize the pinning scheme according to the system topology 
00]. In controlling the dynamics of coupled map lattices, it 
is found that, with the same pinning effort, the controllabil- 
ity can be significantly improved if the pinnings are arranged 
in an asymmetric fashion |g]. The underlying mechanism for 
this improvement is that, by asymmetric pinning scheme, the 
topology of the lattice network is actually modified from sym- 
metric to asymmetric, which causes a more uniform suppres- 
sion of the unstable modes |6]. Still in lattice network, based 
on the special property of coupling gradient, it is even found 
that the spatiotemporal chaos can be successfully controlled 
by pinning a single node floTl , In this application, to iden- 
tify the location of the pinned node, we need to analyze the 
topology of the weighted and directed lattice 1 10]. 

Comparing to the regular networks, optimization of pinning 
scheme according to topology properties is more significant in 
complex networks fill - mill . For practical networks, a general 
feature is that the number of links (degree) associated to the 
network node has a large variation, e.g. the type of scale- 
free networks (SFN) B17ll . In networks like SFN, the ma- 
jority of the network nodes have small degree, but a few of 
the network nodes could have distinctively large degree. No- 
ticed the important roles of the large-degree nodes played in 
other types of network dynamics, e.g. network synchroniza- 
tion ni8fl . it is natural to expect that the spatiotemporal chaos 
of complex networks will be better controlled by pinning the 
large-degree nodes than pinning the small-degree ones or pin- 
ning randomly. To check this, in Ref. ill 111 the authors com- 
pared the performance of two pinning schemes on SFN: large- 



degree pinning (LDP) and random pinning. In LDP, the nodes 
are pinned by the decreasing order of their degree, while in 
random pinning the nodes are randomly selected. The authors 
found that, with the same number of pinned nodes, the LDP 
strategy has a clear advantage over the random pinning strat- 
egy. After the work of Ref. HHl . pinning control of complex 
networks has received continuous interest, and the advantage 
of LDP has been well addressed lll3l4l6ll . 

In a practical situation, due to the limited pinning strength, 
to control a large-size complex network, a number of the net- 
work nodes should be pinned simultaneously ifTTI IT2I1 . For 
instance, with a reasonable pinning strength (even with the 
order of the dynamics variables), it is shown that the number 
of pinned nodes necessary for controlling a SFN to a homo- 
geneous state could take up to 30 or 50 percent of the total 
network nodes ill ill . With such a significant fraction of nodes 
be pinned, it is questionable whether the LDP strategy will 
be still the best option, as now many of the pinned nodes do 
not have large degree. In this paper, we will re-evaluate the 
performance of LDP by comparing it with its opposite case, 
namely the small-degree pinning (SDP) strategy. In contrast 
to LDP, in SDP the network nodes are pinned by the increas- 
ing order of their degree. Very interestingly, we found that 
when the fraction of pinned nodes is large enough (depending 
on the network details, this fraction could be ranging from 30 
to 50 percent), the SDP strategy, which has been tacitly be- 
lieved as having the worst performance, will outperform the 
LDP strategy. We demonstrate this phenomenon in various 
complex networks, and investigate its underlying mechanism 
by a simplified network model. 

Our model of network pinning is the following 
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[H(x J )-H(x,)]+r ? <WH(x T )-H(x 4 )] 

(1) 

where i,j = 1, 2, . . . , N are the node indices, F is the func- 
tion of the node dynamics, H is the coupling function, e and r\ 
represent the uniform coupling strength and pinning strength, 
respectively. The network structure is captured by the adja- 
cency matrix A = {a^}, in which = 1 if nodes i and j 
are directly connected, and = otherwise. The degree of 
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node i thus reads fej = Ylj=i a ij- The target orbit to which the 
whole network is assumed to be controlled is denoted by xy. 
The total number of pinned nodes is n, and the set of pinned 
nodes is represented by V = {m}. Specifically, if node i is 
pinned in the network, we have <5j m = 1 in Eq. [1], otherwise 
Sim = 0. Without loss of generality, we set the target (con- 
troller) to be having the same dynamics as the network node, 
i.e. xt = F(xx). 

Regarding the controller as an additional node of the net- 
work, we then are able to treat the pinning problem under the 
framework of network synchronization 11911 . In the enlarged 
network, the controller has degree n, and is unidirectionally 
coupled to the pinned nodes with strength r\. So the enlarged 
system can be regarded as a special example of weighted and 
directed network. The enlarged network can be written as 

/ a n aw ■ ■ ■ o>in «>i(jv+i) ^ 

_ an ai2 ... Q>iN w i(N+l) 
IJVl &N2 ■ ■ ■ O-NN Wn(N+1) 

\ ... / 

in which the controller is represented by the node of index 
N + 1. The pinnings are represented in the last column of W, 
where tu^jv+i) = r] if i £ V, and uvjv+i) = otherwise. 
Based on the method of master-stability function (20II21I1 . we 
know that whether the network can be "synchronized" to the 
controller is determined jointly by the node dynamics, F, and 
the eigenvalue spectrum of the coupling matrix G = W — DI. 

Here, D — (di, d 2 , ■ ■ ■ , d N +i) T (di = w ij is tne cou- 

pling intensity of node i) and I is the identity matrix of dimen- 
sion N + 1. Previous studies have shown that, for the general 
types of node dynamics, the network is more synchronizable 
when the spread of the eigenvalue spectrum is narrower 11201 - 
l22tl . In particular, let = Ai > A2 • • ■ ^n+i be the eigen- 
value spectrum of the coupling matrix G. Then the smaller 
the eigenratio R = A7V+1/A2, the more likely synchronous 
dynamics is to occur on the network. Since here we treat the 
pinning problem as a special case of synchronization, smaller 
eigenratio thus also means higher controllability. 

Having specified the pinning model, we now compare the 
performance of the two pinning strategies, namely LDP and 
SDP, based on numerical simulations. Without loss of gen- 
erality, we generate a complex network by the standard BA 
algorithm B17I1 . which consists of N = 1000 nodes and has 
average degree (fc) = 8. The node degree of the generated net- 
work follows roughly a power-law distribution p(k) ~ fc -7 , 
with 7 = 3. The LDP and SDP strategies are implemented 
as follows. We first rearrange the nodes by a descending 
order of their degree, i.e. ki > k 2 ■ ■ ■ > fcjv. Then, in 
LDP the n largest-degree nodes are pinned all together, i.e., 
nodes of index 1 < m < n are pinned; while in SDP, 
the n smallest-degree nodes are pinned i.e., nodes of index 
N — n + 2 < to < N + 1 are pinned. The eigenratio of 
the two pinned networks are denoted by R L (LDP) and R s 
(SDP). The variations of R L and R s as a function of the frac- 
tion of pinned nodes, / = n/N, are plotted in Fig. 1. In- 
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FIG. 1 : (Color online) With pinning strength r\ = 25, the comparison 
of the performance between the LDP and SDP strategies in control- 
ling a scale-free network. Inset: the variations of the largest non-zero 
eigenvalue as a function of the fraction of pinned nodes. Each data is 
averaged over 20 network realizations. 



terestingly, we find that, when the number of pinned nodes 
is large enough, R s is smaller than R L , i.e., the SDP strat- 
egy has a higher performance than the LDP strategy. More 
specifically, there exists a critical fraction, f c w 0.4, in the 
number of pinned nodes. When / is small, as tacitly believed, 
the LDP strategy has the higher performance than the SDP 
strategy. However, as / increases, the advantage of LDP over 
SDP is gradually narrowed and, at the critical value f c , the 
two strategies give the same performance. After that, the SDP 
strategy will outperform the LDP strategy. Moreover, as / in- 
creases from f c , the advantage of SDP over LDP is gradually 
enlarged. The difference between the two strategies reaches 
its maximum at / = 1 — 1/N, where only the small-degree 
(largest-degree) node is un-pinned in LDP (SDP). Of course, 
at the point / = 1 the network is globally controlled, and the 
two strategies give the same performance again. (The switch 
of control performance is also observed in the variation of A2, 
which characterizes the system controllability for the situation 
of non-bounded MSF function lf22ll .) 

Is the above phenomenon general for complex networks? 
To answer this question, we have compared the performance 
of the two strategies on a variety of networks, including 
changing the degree exponent, the average degree, the net- 
work size and the pinning strength of the network model used 
in Fig. 1. In exploring the influence of the degree exponent 7, 
we have adopted the network model of Ref. 112311 . where 7 can 
be increased from 3 to a large value by tuning a parameter. In 
doing this, the network topology is changing from heteroge- 
nous to homogeneous gradually. The variation of the critical 
fraction f c as a function of 7 is plotted in Fig. 2(a), where it is 
found that the value of f c is monotonically decreased. Mean- 
while, with the increase of 7, the performance difference be- 
tween the two strategies, Ai? = R L — R s , is narrowed, e.g, 
the variation of Ai? at n = N — 1 [Fig. 2(b)]. The narrowed 
difference in Fig. 2(b) seems to implies that, comparing to 
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FIG. 2: The other parameters are the same as in Fig. 1. By varying 
the degree exponent, the variations of (a) the critical fraction f c and 
(b) the performance difference between SDP and LDP at n = N — 
1. (b) and (c), the variations of f c and ARn-i as a function of 
the network average degree. Each data is averaged over 20 network 
realizations. 



homogeneous networks, heterogeneous networks are more fa- 
vorable for SDP. To check the influence of the average degree, 
we have increased the value of (fc) from 6 to 80. The vari- 
ations of the critical fraction and the performance difference 
between the two strategies are plotted in Fig. 2(c) and (d), re- 
spectively. It is seen that, as (fc) increases, both f c and AR 
are decreased. Besides degree exponent and average degree, 
we have also checked the influence of the network size, and 
found that, as N increases from 1000 to 5000, the values of 
f c and AR are decreased too (not shown). Noticing the fact 
that either the increase of (fc) or the decrease of N leads to a 
homogeneous network structure, these additional simulations 
therefore confirm the finding of tuning 7, i.e., the SDP strategy 
for more favorable in controlling heterogeneous networks. 

The existence of f c can be mathematically analyzed, based 
on the model of star networks. When not pinned, the star net- 
work has three distinct eigenvalues: Ai = 0, A^v = — N, and 
Aj = —1 for i = 2, . . . , N — 1. When pinned, the eigen- 
value spectrum will be modified, reflected as the decreases of 
Ai, Ajv and a few of other eigenvalues. Specifically, when n 
nodes are pinned, in LDP the number of modified eigenvalues 
is n + 1, while in SDP this number is n + 2. Firstly, let us 
consider the case of n = 1 and compare the performance of 
the two strategies. For LDP, only the central node is pinned, 
and the two modified eigenvalues, Af and A^, are governed 
by equation 

A 2 + (N + r?)A + 77 = 0. (2) 

We thus have Af >N = [-(N + rj) ± ^(N + rj) 2 - 4^/2. For 
SDP, the pinned node could be any of the A — 1 peripheral 
nodes. The three modified eigenvalues, Af , Xf f _ 1 and Ajy, 
are governed by equation 

-(A 2 + AA + 1)?7=(A + A)(A + 1)A. ( 3 ) 



By Cardan formula, we have Af = — (A + 1 + 77 + tuip 1 / 3 + 
cj 2 g 1/3 )/3 and A^ = -{N + 1 + 77 + p 1/3 + g 1/3 )/3, 
where p = (a + Va 2 -46 3 )/2, q = (a - Va 2 -46 3 )/2, 
a = 2(A + 1 + rj) 3 - 9A(A + 1 + r?)(l + 77) + 27rj, 
b = (N + 1 + T]) 2 - 3A(1 + rj), and Wi , 2 = -0.5 ± 0.5%/3i 
Previous studies have shown that, to control the dynamics of 
a complex network successfully, the pinning strength should 
be sufficiently large, i.e., rj >> 1 llll1 - [l3ll . Under this con- 
dition, the above eigenvalues are approximated as Af ~ —1, 
~ -(V + N), Af w —1/N, Xff w -(77 + I). Apparently, 

R s = Xff/Xf > R L - Aft/Af, (4) 

i.e., LDP has the higher performance than SDP when 1 node 
is pinned. 

We next compare the performance of the two strategies for 
the situation of n = N — 1, i.e. at the end of the comparison. 
For LDP, this means that except one peripheral node, all other 
nodes on the network are pinned. In this case, the eigenvalues 
Xj^_ 1 and Xj^_ 1 are governed by equation 

- [(JV + 2A)?7 + r; 2 ](A + 1) + V = (A + N)(X + 1)A. (5) 

While for SDP, only the central node is un-pinned. In this 
case, the eigenvalues Xf f _ 1 and A^)^ are governed by equa- 
tion 

A 2 + (N + r])X + Nr) = 0. (6) 

Under the condition 77 >> 1, these eigenvalues are approx- 
imated as Af w -1, A^ w -(77 + 1), Af w —N, X% « 
— (77 + 1). So, when — 1 nodes are pinned, we have 

R s = X%/Xf <R L = A^/A[, (7) 

i.e., SDP has the higher performance than LDP. Since the sit- 
uations n = 1 and n = N — 1 stand as the two ends of the 
comparison [Fig. 1], Eqs. [4] and [7] thus guarantee the exis- 
tence of at least one crossing in the performance curves. 

The above phenomenon of performance switching can be 
physically explained as follows. Still take the star-network 
as the model. When only one node is pinned, by pinning 
the largest-degree node, the pinning signal can be efficiently 
propagated to the peripheral nodes, as the hub-node has the 
shortest average-node-distance. For this reason, LDP will 
have a higher performance than SDP (see Eq. [4]). How- 
ever, this picture is changed when N — 1 nodes are pinned. 
Under a strong pinning strength (77 >> 1), in LDP the N — 1 
pinned nodes have been well controlled to the target, and they 
will influence the un-pinned peripheral node via a coupling of 
strength e; on the other hand, in SDP all the peripheral nodes 
have been controlled, and they influence the central node to- 
gether via a joint coupling of strength (N — l)e. Therefore 
in SDP the coupling between the pinned part (the peripheral 
nodes) and the un-pinned part is in fact amplified by N — 1 
times, therefore making SDP superior to LDP. Having under- 
stood this mechanism, we are able to predict that, as the size of 
the star network increases, the advantage that SDP over LDP 
will be enlarged. Actually, this point can be also drawn from 
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FIG. 3: (Color online) Comparison of performance the two pinning 
strategies, LDP and SDP, for (a) the power-grid network of the west- 
ern United States and (b) the Internet at the Autonomous level. 



Eq. [7], which implies that AR = R L - R s oc N. As the de- 
gree of heterogeneity of star networks can be characterized by 
the network size, this relation thus implies from another angle 
the superiority of SDP in controlling heterogeneous networks. 
It is worth noting that a similar phenomenon about the role of 
small-degree nodes has been also observed in network syn- 
chronization, where it is found that, to achieve global network 
synchronization, it is usually the smallest-degree nodes that 
are more important I 24I1 . 

How about the realistic networks? The specific examples 
we have checked are the electrical power grid of the western 
United States 12511 and the Internet at the autonomous level 
l26ll . The power-grid network consists of N — 4941 nodes 
and has average degree (k) w 2.67. The variations of R as 
a function of n for both pinning strategies are plotted in Fig. 
3(a). Interestingly, it is found that, different to the results of 
the ideal models [Fig. 1], the value of R is decreased in a 
non-smooth fashion and the two curves are crossed at several 
places. We attribute this phenomenon to the inherent commu- 
nity structures in the power-grid network (a detail analysis on 
the non-smooth variation of R will be presented elsewhere). 
The Internet we have employed consists of N = 6474 nodes 



and having an average degree (k) w 3.88. The variations of 
R as a function of n for both pinning strategies are plotted in 
Fig. 3(b). Again, we find a non-smooth decrease of R. For 
the Internet network, we find only one crossing between the 
two curves, which is located at about n = 2236. 

The finding that heterogenous networks can be better con- 
trolled by pinning the small-degree nodes might have some 
potential applications, say, for example, in synchronizing sen- 
sor networks or controlling the consensus of multi-agent sys- 
tems. In a sensor network, due to terrain complexity the sen- 
sors are distributed in a non-uniform fashion: most of the sen- 
sors have very few neighbors, but a few sensors could have 
many neighbors 12711 . That is, the sensor network has a het- 
erogenous degree distribution. To synchronize the timers of 
the sensors, a general approach is to assemble a fraction of 
the sensors with a signal receiver, which receives signal from 
the satellite and reset the timer of the sensors accordingly. For 
those sensors without the receiver, their timers are updated by 
their neighbors. Here an important question is that, for the 
given number of sensors having signal receiver, how to place 
them properly in order of a higher performance for timer syn- 
chronization? The study of the current paper suggests that, if 
the number of receiver-assembled sensors are large enough, it 
may be better to place them in the areas of sparse sensors. An- 
other application could be the control of consensus in multi- 
agent systems l28ll . To make the mobile agents move in a 
synchronous fashion, in engineering one practical approach is 
to introduce some "leading" agents into the system. Differ- 
ent to the normal agents, these leading agents have a prede- 
fined motion, but they will influence the motion of the normal 
agents via couplings (behave like the pinning controller). Our 
study suggests that, to control the multi-agent system more ef- 
ficiently, we may choose to put the "leading" agents into the 
small swarms, instead of the tacitly believed large swarms. 
(From the network point of view, large and small swarms 
could be regarded as possessing large and small degree, re- 
spectively.) 

In summary, we have compared the performance of two op- 
posite pinning strategies in controlling complex networks. It 
is found that, when a significant fraction of the network nodes 
are pinned, the pinning of the small-degree nodes could have 
a higher performance than the large-degree ones. We have ex- 
plained this phenomenon by a model of simplified network, 
and discussed the dependence of this phenomenon to the net- 
work parameters in detail. These findings, hopefully, will pro- 
vide a new angle to the pinning control of complex networks, 
and apply to practical situations. 
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